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Abstract 


The  primary  result  of  this  research  effort  has  been  the 
development  of  an  image  algebra  that  can  serve  as  the  foundation 
of  a  common  algebraically  based  image  processing  language,  in 
comparison  to  other  existing  image  algebras,  this  algebra  is 
capable  of  expressing  common  image  processing  algorithms  and 
transforms  in  terms  of  its  operators. 
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Introduction 


Current  image  processing  algorithm  development  is  not  based 
on  an  efficient  mathematical  structure  that  is  designed 
specifically  for  image  manipulation,  feature  extraction  and 
analysis.  In  general,  each  researcher  develops  his  own  set  of 
ad-hoc  image  processing  tools,  thereby  increasing  research  and 
development  costs  accordingly.  The  vast  increases  in  image 
processing  activities  throughout  the  military,  industrial  and 
academic  communities  are  resulting  in  an  immense  proliferation  of 
different  operations  and  architectures  that  all  too  often  perform 
similar  tasks.  There  are  probably  as  many  image  processing 
languages  as  there  are  architectures,  and  all  of  them  differing 
in  capabilities. 

In  view  of  this  ever-increasing  diversity  of  image 
processing  architectures  and  languages,  the  principal 
investigator  proposed  to  develop  a  standard  image  processing 
algebra  to  serve  as  a  mathematical  basis  for  a  common  image 
processing  language.  The  relational  formalism  of  an  algebraic 
image  processing  language  would  constitute  an  invaluable  aid  in 
the  design,  development,  optimization  and  testing  of  image 
processing  algorithms  and  hardware  configurations. 

This  research  effort  succeeded  in  defining  a  universal  ir, age 
algebra  that  could  serve  as  the  origin  from  which  a  common  image 
processing  language  could  evolve.  As  compared  to  other  existing 
image  algebras,  no  problems  ha>e  been  encountered  in  translating 
common  image  processing  opera.  into  the  language  of  this 
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algebra.  Because  of  this  enormous  success,  research  efforts  are 
now  continuing  on  four  fronts:  (1)  Further  theoretical 
development  of  the  algebra;  (2)  Formulation  and  compilation  of  a 
new  image  processing  language  based  on  this  algebra;  (3) 
Optimization  of  FLIR  algorithms;  (4)  Design  of  reconf igurable 
VLSI  architectures  for  image  processing  based  on  this  algebra. 


Technical  Aspects  of  the  Research  and  Research  Results 

This  research,  under  Air  Force  grant  No.  AFOSR-83-0065 ,  has 
been  involved  with  the  development  of  a  rigorous  formulation  of  a 
mathematical  foundation  for  image  processing  algorithms  and 
operations.  The  research  proceeded  along  the  following 
guidelines; 


(1)  Investigation  of  existing  algebraic  structures  for  image 
processing . 

(2)  The  cataloguing,  according  to  task  similarity  of 
existing  image  processing  operations. 

(3)  Investigation  of  the  relationships  between  the  basic 
components  oI  the  catalogued  operations. 

(4)  Extraction  of  a  minimal  set  of  operators  in  order  to 
form  an  algebra  capable  of  expressing  all  image 
transform  operations. 

(5)  Establishment  of  some  basic  relationships  and  theorems 
governing  the  algebra. 
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(6)  Investigation  of  the  algebra's  potential  to  serve  as  a 
foundation  of  a  common  algebraically  based  image 
processing  language. 

The  investigation  of  existing  image  algebras  turned  out  to 
be  somewhat  disappointing.  Extracting  the  mathematical  and 
conceptual  core  of  existing  algebras  (ref.  6,11,12),  resulted  in 
only  one  structure  mathematicians  would  dare  call  an  "algebra". 
This  structure  is  equivalent  to  the  Minkowski  algebra  of  sets 
(4),  Although  the  literature  abounds  with  so  called  new 
techniques  -  i.e.  erosions,  dilation  closings,  openings,  rolling 
ball  algorithms,  etc.  -  the  algebraic  relationship  provided  by 
these  techniques  can  all  be  found  or  easily  derived  from  the 
algebraic  relationships  in  (4). 

Even  though  many  neighborhood  operations  can  be  expressed  in 
terms  of  the  Minkowski  algebra,  the  algebra  is  extremely  limited 
in  performance  (6),  incapable  of  expressing  global  transforms  and 
various  neighborhood  transforms  and  can,  therefore,  never  serve 
as  a  universal  image  algebra.  It  became  clear  at  the  initial 
stage  of  this  research  that  a  new  algebra  needed  to  be  defined. 

In  order  to  accomplish  this  task,  cataloguing  and  investigating 
existing  image  operations  became  a  necessity. 

Because  of  the  "one  year"  tim-  constraints  some  image 
transforms  were  more  thoroughly  investigated  than  others.  In 
particular,  emphasis  was  placed  on  linear  transforms,  non-linear 
smt  ithing  and  enhancement  techniques,  edge  detection  schemes, 

1C  sf>gmentat ion  and  background  removal.  W.K.  Pratt's  book  on 
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digital  image  processing  (7)  proved  to  be  the  most  valuable 
resource  for  this  particular  task. 

Investigation  of  the  basic  components  of  the  various  image 
transforms  made  it  clear  that  we  were  dealing  with  only  four 
basic  operations,  namely  two  arithmetic  and  two  logic  operations. 
However,  these  four  operations  are  applied  differently  in 
cellular  (neighborhood)  processing  than  in  non-cellular 
processing,  resulting  in  an  algebra  of  eight  binary  operators. 

The  definition  of  these  operators  can  be  found  in  Appendix  2  and 
details  concerning  the  algebra  are  given  in  (8). 

In  order  to  demonstrate  the  potential  of  this  algebra  to 
serve  as  a  basis  for  a  common  image  processing  language,  we 
showed  that  the  algebra  is  capable  of  describing  all  commonly 
used  image  processing  functions.  We  considered  such  diverse 
processing  operations  as  linear  transforms  (Fourier  transforms, 
Walsh  transforms,  etc.),  non-linear  filtering  and  enhancement 
techniques,  a  variety  of  well-known  edge  detection  algorithms, 
gray  scale  averaging,  thresholding,  and  histogram  equalization. 

In  order  to  minimize  the  number  of  operations  in  the  algebraic 
formulation  of  algorithms,  we  also  established  some  basic 
relationships  governing  the  algebra.  As  mentioned  earlier,  no 
great  difficulties  were  encountered  when  translating  image 
operations  into  our  algebraic  formulation,  A  special  bonus  of 
the  translation  task  was  the  discovery  of  some  new  and  powerful 
image  enhancement  techniques  (9). 
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Discussions  with  my  colleague  Dr.  S.  Chen  of  the  National 
Science  Foundation  led  to  the  discovery  that  the  algebra  is  an 
"image  processing  machine"  in  the  abstract  sense,  and  can  thus  be 
used  to  define  and  model  real  architectures.  Building  on  this 
idea,  we  defined  a  language  based  systems  architecture  where  the 
algebraic  algorithms  are  expressed  as  data  flow  graphs  that  are 
mapped  to  a  reconf igurable  distributed  system  (3).  In  view  of 
recent  advances  in  VLSI  technology,  such  architectures  are  now 
feasible. 

In  our  system,  the  user  inputs  image  data  through  a  front- 
end  computer  to  a  distributed  network  which  leads  to  various 
operation  modules.  The  active  operation  modules  drive  parallel 
processing  elements  that  carry  out  the  elementary  algebraic 
operations  and  transformations.  Configurations  for  our  variable 
neighborhood  definition  are  formed  through  the  control  of 
arbitration  networks.  Modularity  and  redundancy  will  enable  the 
system  to  be  fault  tolerant  and  expandable. 

The  main  improvement  of  this  system  over  some  existing 
architectures,  such  as  the  cy tocoraputers  or  the  CLIP  series,  is 
the  arlity  to  handle  variable  neighborhoods  and  to  perform 
certa:*’  image  processing  algorithms  in  parallel  that  are  not 
feasibi*^  on  current  cellular  array  computers. 

P  .ide  benefit  of  this  research  was  the  development  of 
software  for  the  VAX-11/780  computer  to  enable  the  printing  of 
gray  level  images  on  standard  dot-matrix  printers.  The  extremely 
low  resolution,  spatial  distortion  and  slowness  of  the  print 
routines  provided  by  Government  furnished  print  routines 
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necessitated  this  development.  This  software  will  be  made 
available  to  AFATL  personnel  at  Eglin  AFB.  Appendix  1  provides 
examples  of  the  improved  image  displays. 

Finally,  this  one-year  research  effort  resulted  in  four 
publications  (1,2,9,10)  and  four  invited  talks  and  lectures.  In 
October  1983,  the  principal  investigator  was  an  invited  speaker 
and  session  chairman  at  the  IEEE  International  Conference  on 
Computer  Design  held  in  New  York.  At  the  1984  annual  Spring 
meeting  of  the  Mathematical  Association  of  America's  Florida 
section  in  Tampa,  an  invited  lecture  was  given  on  the  connection 

between  digital  topology  and  the  image  algebra.  Two  talks 

« 

concerning  the  algebra  were  given  in  April  1984,  one  at  the 
Conference  on  Intelligent  Systems  and  Machines  in  Rochester, 
Michigan,  and  the  other  at  the  1984  Southeast  Regional  ACM 
Conference  in  Atlanta,  Georgia. 


/ 


Summary  and  Recommendations 

We  constructed  an  algebra  for  image  manipulations  consisting 
of  eight  binary  operators.  All  image  processing  techniques 
investigated  during  this  research  were  expressible  in  terms  of 
this  algebra,  and  it  is  our  opinion  that  most,  if  not  all, 
current  image  manipulation  techniques  are  translatable  into  the 
language  of  this  algebra.  However,  this  needs  to  be  more 
thoroughly  documented. 

Since  very  little  is  known  about  algebraic  structures 
containing  more  than  two  operators,  theorems,  corollaries, 
identities,  and  laws  concerning  compositions  of  different 
operations  need  to  be  established.  Such  laws  will  not  only  be 
useful  in  terms  of  algorithm  simplification  but  also  provide 
deeper  insight  and  a  better  understanding  of  the  algebra. 
Following  the  establishment  of  such  relationships  and  laws,  a 
natural  "next  step"  could  be  the  application  of  these  algebraic 
relationships  for  the  optimization  and  testing  of  Government 
furnished  FLIR  algorithms,  and  also  as  an  aid  in  the  development 
of  new  image  processing  algorithms.  The  importance  of  such  a 
fully  developed  algebra  with  regard  to  military  applications 
cannot  be  overestimated. 
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AESTEACr 

Iq  this  paper  we  define  an  leugc  algebra 
whose  operators  serve  as  the  basis  of  a  new  Inage 
processing  language*  In  comparison  to  other 
existing  Image  algebras,  this  algebra  Is  capable 
of  cspresslng  nost  coamoo  Icage  processing  algor* 
ithas  and  cransforaa  In  ceres  of  Its  operators* 
The  development  of  this  algebra  has  been  Influ* 
enced  by  Che  archlceeturea  of  oasslvely  parallel 
luge  processlog  systeas* 


1*  IKTBODUCTIOH 

About  25  years  ago,  Unger  (16)  proposed  that 
algorithm  for  luge  processing  and  analysis  could 
be  lepleaentcd  In  parallel  using  "cellular  array* 
cespucers*  Recent  advances  In  VLSI  technology  now 
permits  the  teallzstioo  of  such  array  computers* 

A  detailed  destrlpcloo  concerning  the  diversity 
and  genealogy  of  cellular  array  computers  can  be 
found  lo  (10}« 

for  our  purposes  It  suffices  to  observe  that 
cellular  architectures  Icpleoent  variations  of  von 
Seuaann's  autooacon  (15)* 

S.ASA's  massively  parallel  processor  or  KPP 
(2),  and  Che  CLIP  series  of  coaputers  developed  by 
Duff  (A)  represent  Che  classic  eabodlsenc  of  von 
Neumann's  original  autoeaton.  The  CLIP^  which 
constitutes  Che  latest  In  the  series  of  CLIP 
coopuCera,  consists  of  an  array  of  9216  (96x96) 
processors.  laploySng  VSLI  technology!  sets  of 
eight  processing  elements  are  Integrated  on  a 
alogic  chip. 

The  KPP  also  lncr.r<ttes  eight  processing 
elements  per  chip  in  .  .nsseoblage  of  128x132 
processing  elements.  i  r.  distinction  to  the  CLIP, 
where  each  processing:  .  .  •  “»>nt  has  the  capability 
of  CocsuQicatiog  wit:  eight  Icmedlale  neigh* 

bors,  an  KPP  processl:.  liient  has  connection  to 
only  four  immediate  r.«:.':H;rs  as  Indicated  by  the 
solid  lines  In  Figure  1. 
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trllulir  i»*9»  procrssir.;  iulei*ton  of  idrntici)  pr6C2is9rt 
kith  numl  rin^hbor  ct^^rction 


Figure  1 


Using  these  types  of  hardwired  coemunlcatloD 
links  between  neighboring  processors,  each  pro* 
tesBor  la  responsible  for  one  pixel  Cor  one  ele* 
cent  of  the  Image),  and  la  capable  of  performing 
opersclons  on  the  image  via  its  connur.l  cat  Iona 
links.  These  local  operations  can  be  expressed  in 
terms  of  neighborhood  operators  or  i-orhood 

funtlons  and  are  perlorsed  In  paralle.  u:.  the 
whole  Images  and  nelghborbooda  (l.e.  f.i -Images 
Induced  by  the  local  windows).  In  this  sense  then 
array  processors  impose  a  natural  alg^rra  on  the 
set  of  Images  and  window  configurations. 

Several  Image  algebras  employinp  ::.fse 
concepts  already  exist.  Among  these,  .  .:e  are 
only  three  Chat  maChesatlcians  would  caH 
'algebras"  (7,13,1A).  However,  desplti  thclr 
profound  accocpl  i  shment  B,  these  algctrar.  are  not 
capable  of  expressing  most  comaor  im^ft  processing 
operations  such  as  Fourier  transformations,  gray 
scale  averaging,  and  various  edge  Cit^ction 
techniques.  In  fact,  the  failure  of  these 
algebras  to  express  a  fairly  straight  forward  U.S* 
Government  furnished  FLIR  algorithm  has  tten  well 
documented  (6). 
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to  contr««tt  Image  algebra  developed  by 
(hit  author  la  capable  of  eapreaalng  moat  cobaoq 
iaagc  processlog  operations  In  tenaa  of  Ita 
operators*  The  devclopacDt  was  aotlvatcd  by  the 
Air  Force*!  need  for  translating  image  processing 
algorithms  Into  a  cosmoo  mathematical  language  for 
performance  characterization*  docusentaclon*  and 
algorithm  simplification* 

In  the  oest  section  we  provide  a  rigorous 
mathematical  deflnltloo  of  this  algebra*  an* 
dowlng  It  with  sufficient  flexlbllty  for 
implementation  on  future  reconf igurable 
neighborhood  computers  (3)  as  well  as  com 
^eocional  serial  image  processing  machines* 


2*  rUKDAKEKIAL  TU>1S 

Ueaceforth*  Z  and  C  shall  denote  the  sets  of 
integers  and  complex  numbers*  respectively* 
Although  we  could  Just  as  well  have  used  the  aet 
of  reals  instead  of  complex  numbers*  we  obtain  a 
marhcBatlcally  more  useful  and  extensive  structure 
by  employing  the  latter* 

DCyiNlTION  2.1. 

(1)  S  -  {(x.y.z):  x.yCZ,  xec). 

(2)  P(S)  -  (A:  AC  s}. 

The  power  set  P(S)  will  represent  our 
universe  oi  discourse*  In  particular*  Images  will 
be  viewed  as  elements  of  P(S)* 

Henceforth  s  -  <x,y,z)  and  s'  •  <x*,y',x*) 
will  denote  elements  of  S  and  A*  8  subsets  of  S* 

DEnKTlOff  2.2* 


(1)  s  and  s*  are  said  to  be  related,  denoted 
by  •  —  sf,  if  X  •  X*  and  y  •  y**  otherwise  s  is 
not  related  to  s*,  which  Is  denoted  by 

a  s*. 

(2)  Aj  -  {a  €  A:  a  *  b  for  any  bcfi) 

and  A^*  [a  c  A:  a  ^  b  for  sot>;  t  «■'  B }« 

(3)  A  la  related  to  B*  denoted  by 
A  -  B  If  A®  -  A  xnd  B*  -  B. 


(^)  A  called  an  1  m.i If  whenever  a,b€  A 
aod  a  b,  then  a  •  b. 


If  c  la  a  complex  number  then  the  magnitude 
of  c  will  be  denoted  by  |c|  and  die  real  part  of  c 
by  r(c).  Given  two  complex  nonbers  c  and  c**  we 
define  the  maxi mum  and  el ni mum  of  c  and  c* 


as  c7c*  -  max{r(c),r(c')) 

and  cAc*  •  oln  |r(c) * r(c* ) ) *  respectively* 

DEFINITION  2*3 

(1)  The  gray  level  sum  of  s  and  •*  is  defined 
as  s(+)i*  •  (xTyT*^*  /  tbs 

gray  level  product  as  b(x)s  •  {x*y*sc*}* 

(2)  The  ttsxlmutn  of  s  and  s'  is  defined  as 
s(V)s*  ■  (x,y*zVz*)  and  the  the  minimum  as 
s(A)s’  -  (x*y,zAz*>. 

(3)  If  f  is  a  real  or  complex  valued  function 
on  C,  Chen  the  pixel  functionj^  Induced  by  f  la 
defined  as  _f^(a)  *■  (x,y*f(z)).  In  particular* 
we  denote  the  absolute  value  or  magnitude  of  a 

by  ja|  •  (x«y*|z|),  and  exponentiation  and 
scalar  tajl t Ipllcatlon  by  a  complex  number  c  by 
•  (**y**^)  ci  -  (x*y*cz)*  reapectlvely* 
For  finite  aubseca  of  S  we  alao  define 


the 

follDulng  four  operations 

(4) 

*(+)A  - 

«(♦) 

f.tA>‘ 

(5) 

s(x)A  - 

s(x) 

^!a>* 

(6) 

x(V)A  - 
> 

s(7)( 

.fA>* 

(7) 

«CA)a  • 

a(A)( 

Here  the  operation  (**-)#  means  to  adi^  using  (-^)* 
all  of  a  €  A. 

Several  comments  are  in  order*  First* 
observe  that  the  operations  defined  In  (1)  and  (2) 
are  not  cossutacive*  To  further  clarify  (3)* 
consider  examples  such  as  stn(s)  •  (x»y*slnz) 
exp(s)  -  (x,y,exp(z))  and  Ln(s)  •  (x*y,lnz)* 
Finally*  note  the  order  in  which  the  a's  are  added 
Id  (A)  Is  Inssaterlal  since  s  Is  added  on  the 
"left"*  Thus*  if  A  -  {a*b*c)*  then  s(‘*^)A  - 
B(+)a(+) b(+)  c  •  s(+)c<+)aC  +  )b.  The  fvame 
bbservatlon  holds  for  operations  (5)  through  (7)* 

DEFINITION  2.* 

A  nc ; rhborhood  f unct Ion  or  nelrhborhood  for  S 
is  a  function  h:S  •*  P(S)a  The  mathematical  Image 
N(s)  of  a  point  e  e  S  is  called  the  N  neighborhood 
of  6  or,  simply*  a  neighborhood  of  s.  The 
restriction  of  N  to  a  subset  A  of  S  will  be 
denoted  by  and  Is  called  a  neighborhood  for 
A.  The  deleted  r.clr*  berhood  K*(s)  Is  defined  aa 
N*(b)  •  Ms)  -  \s 

We  are  now  Ir  the  position  of  defining  a 
universal  algebra  whicli  operates  on  subsets  of  S* 


3.  IMAGE  ALCEBBA 

As  earlier*  the  operands  of  our 
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Tbe  algebr*  C»  defined  in  the  previous 
•rction»  Initially  evolved  fro*  Che  (our  neighbor** 
hood  operetlone  (Definition  3.3)  at  a  need  (or 
aathexat leal ly  dcacrlblng  liuge  proceatlng  rout* 
Inee  that  arc  ^natural*  to  cellular  archlcecturcs* 
Addicloo  of  the  rcoalolng  four  arlthaecic-loglc 
operatort,  however,  ylclda  a  aore  flexible 
algebra*  la  fact,  the  extended  algebra  provldca  a 
unlfora  aethod  (or  deeerlblng  aoic  atandard  Iraec 
traneforaa  and  laage  proceatlng  techniques  in 
teraa  of  algebraic  foroulae  Involving  only  imaget 
aad  the  opetaciont  defining  the  algebra*  Thla  la 
accoepllshed  by  analyzing  the  basic  coaponenta  and 
operacloos  constlcutlog  a  given  aanlpuletion  or 
transfora  aod  then  translating  then  into  the 
language  C.  Due  to  Halted  space,  we  precenc 
traoslatiooa  of  only  a  few  well-known  technlquea. 
oaltcing  Bore  coopllcated  algorithms,  proofs  and 
vcrlf Icatiooa.  all  of  which  can  be  found  on  (»#). 


4.1  Eoce  DKTECriOII 


The  edge  decectloc  Shi  eJ'-hai.ceaent  technique* 
described  In  (8)  ere  easily  trenslsced  Into  the 
language  of  u.  Ue  provide  two  standard  exaoples. 


Defining  the  deleted  neighborhood  function 


by  B‘(a) 


"l+3 . 1*7' 

permits  us  to  express  the  Klrtch  edge  detection 
algorithm  (6)  ms 


K  -  iv{  ^  )|  a(+)b;  *  a|), 

1-0  * 

where  C  deootea  the  enhanced  laage  obtained  froa 
A.  The  logarithmic  edge  detection  scheme  aa 
defined  by  Wallis  (6.pe4B9)  translates  Intoi 


W  -  (l/4)Lnl\*(A^x)(  X  )(M’  )■*)}, 
1-0 


where  A  and  V  denote  t;e  Input  and  output  litages. 
respect Ively. 


4.2  T'iSfSSOLDLNC 


prndently  developed  by  P.E.  Hiller  (7),  J.  Serrs 
(13)  and  S.  Sternberg  (14).  The  two  fundaaenlal 
operatora  of  thla  algebra  correspond  CO  ths 
Hlnkowakl  addition  and  aubtractlon  of  stta  In 
Euclidean  spact  (12)  and  (5).  In  the  Isiags 
procesalng  literature  the  Hlnkowakl  opcratlona  ars 
coaaonly  referrad  to  aa  the  espanalon  or  dilation 
and  Che  erosion  or  shrlnklns  operttors.  It  turna 
out  that  If  N  la  any  neighborhood  configuration 
and  A  .  .■>  Inage,  then  the  dilation  of  A  by  N 
tranlates  Into  A(T)N  and  Che  erosion  of  A  by  N 

Into  A(A)N.  For  exanple.  If  N(a)  •  Hq24a(*) 

and  A  Is  aa  shown  In  Figure  2(a),  then  A(V)N 
and  A(A)N  arc  aa  ahown  In  Figure  2(b)  and  (c)« 

It  follow!  that  our  algebra  generalizes  ths 
Minkowski  image  algebra. 

Two  of  Che  most  basic  and  far  reaching 
comblnaclona  of  the  Minkowski  operations  have 
become  known  as  the  closing  and  opening 
operations.  A  closing  is  an  expansion  followed  by 
a  shrinking  while  an  opening  la  a  shrinking 
followed  by  an  expand.  For  A  and  N  as  In  the  last 
example,  the  doting  of  A  by  N  corretpnds  to 
(A(V)K)(A)N  and  it  shown  In  Figure  2(d).  Observe 
chat  the  result  of  a  closing  la  a  smoother  image, 
with  Che  Interior  'holes*  (zeros)  removed.  In 
many  cases  the  closing  filter  exceeda  the  local 
averaging  filter  on  performance  (1). 


00000000000000 
^0000010000000 
^OOOOIlIlOllOO 
0011 11 11 111000 
lOOOOl  100011000 
.00000111111100 
000000000 i 0000 
00000000000000 

(a) 


000000 1 0000000 
00000111101100 
00111111111110 
oiiiiiniuloo 
00111111111100 
00001111111110 
00000111111100 
OOClOOOOOO  10000 

(h) 


00000000000000 

00000000000000 

00000010000000 

00000100000000 

ODOOOOOOOOOOOO 

00000000010000 

ODOOOOOOOOOOOO 

ODOOOOOOOOOOOO 

(c) 


00000000000000 
00000010000000 
00000111101100 
00111111111000 
00001  nil  11000 
0000011 1111100 
00000000010000 
00000000000000 

(d) 


Figure 

4.5  THK  FOUSXEk  TkAhSFOEM 


Since  (rI*)iA  •  Wa.b.rc):  (a.b.c)  a1,  the 

luge  obtained  from  thresholding  A  at  given  by 

B  -  (rl*).:A. 


A6  a  final  exa-pie  we  examine  the  discrete 
Fourier  transorm  (DFT)  cf  a  nxn  image  array.  The 
usual  definition  of  the  DFT  is  given  by 


4.3  AVXIACDC  niTU 

Syppoae  .V(a)  deaoces  the  window  with  with 
center  pixel  e  used  for  avcr:*qlng  A.  Then 
slace  A(+)N*  *  {a('*')K '  (a):  atA],  the  averaging 

filter  can  be  translated  aa  B  -  n  Na(4)N*),  where 
o  represents  the  number  of  pixels  In  the  window  N 
and  B  the  output  obtained  from  the  Image  A. 

4.4  CEQKtTXiC  flLTEBS 

As  mentioned  In  the  Introduction,  a  fairly 
successful  Boolean  luge  neighborhood  algebra, 
based  on  two  fundamental  operators,  was  Indc- 


n-1  D“1 

F(u,v)  -  Cl/n)  I  I  lf(x,y) 
x-0  y-0 

f>p]f-lT./n)(u,y+vy))  ), 

where  f(x.y)  represents  the  gray  level  at  position 
(x.y). 

In  order  to  express  this  transformation  In 
terms  of  whole  iasges.  we  need  to  define  some 
special  luges  and  neighhorhooda: 
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«lg*br«  wll  be  lAege*  end  nelghborhoodi  (tub* 
end  windows)*  Before  defining  ths 
operscors  of  this  slgebrs*  we  have  need  to  define 
eertslo  specisl  loegea* 

Dt/limON  3.U 

(1)  The  rero  lasge  la  defined 

*•  0  -  {(z.7.t>  «  S:  *  -  o}. 


■» 

A»A  *  •  1*  •nd  a"  “  AiiAx...xA,  wtier*  n  • 
positive  integer  and  the  product  conalats  of 
n^factors*  further  properties  and  theoreos 
concerning  these  operations  can  be  found  In  (10)* 

The  next  aet  of  operations  we  arc  about  to 
define  is  a  special  set  of  neighborhood 
operat lens >  However,  we  point  out  that  the  neigh* 
borhood  N  could  be  replaced  by  any  aubset  B  of  S. 

DmMITlON  3.3 


(2)  The  unit  ioage  la  defined  aa 

I  •  € S!  X  -  i). 

(3)  Z/  f  la  a  real  or  cosplex  valued 
function  oo  C,  then  the  induced 

iaage  function  F  oo  P(S)  Is  defined  as 

F(A)  *  {f^(a);  aCA],  for  each  subset 
A  of  S. 

Here  denotes  the  pixel  function  Induced  by  f* 
Thus  Exp(A)  *  fexp(a) :  a  a|  and  Ln(A)  * 

{ln(a):  a«A}«  In  particular^  the  pagnltude  of 

A  la  given  by  [a|  *  {|*|*  cA  * 

[ca:  a  <£  A I  and  A^  *  (a^:  a  c  A )  *  where  c  €  C* 

Thus,  -A  •  {-a:  a«A}  -  ((K.y.“*)! 

(x.y,z)cA]  jnd  a“*  V  {*"*  :  aCA)  - 
{(x.y.l/z):  (x.y.sllA)*  Note  that  the  latter  la 
defined  only  If  A  •  g* 

Ve  are  oow  ready  to  define  the  first  four 
binary  operators  of  our  algebra* 


Let  N  be  a  finite  neighborhood  of  A*  The 
neighborhood  suo,  raaxloutt  and  alnlnua  of  A  and  H 
are  respectively  defined  as; 


(1)  A(+)N 

(2)  A(x)N 

(3)  A(V)K 

(4)  A(A)N 


{a(+)NCa):  a  e  a) 
(aCx)N(a):  a  C  a) 
{a( V)N(a):  a  c  a] 
{a(A)N(a>;  aC  a). 


The  unlveraal  image  algebra  Is  now  defined  as 
the  pair  Q  -  {J,T)  where  T  •  |+,x,VpJ^(+),(x), 
(?),(A)}  and  J  denotes  the  set  of  Images  and 
nelghborhooda*  Varioua  properties  of  this  algebra 
have  been  explored  In  (10)  and  (11)*  fhe  next 
section  provides  but  a  small  glance  at  the 
potential  of  this  algebra*  The  examples  we  give 
should  also  provide  sufficient  insight  Into  Che 
natural  interaction  of  this  algebra  with  the 
architecture  and  operations  of  cellular  array  ^ 
cospucers*  This  is  particulary  evident  whea 
considering  the  last  four  neigborhood  operations** 


4*  APPLICATIONS 


Deraanoif  3*2 

Let  A  and  B  be  Isuiges  and  C  ■  w  B^*.* 

(1)  The  gray  level  sum  of  A  and  B  Is 

defined  as  A  4  B  *  {a(^)b;  aCA, 
b  €B,  and  a  —  b}  w  E. 

(2)  The  gray  level  product  of  A  and  B 

is  defined  as  AxB  *  {a(x)b:  s£A, 
b  CB*  and  a  **  bj  w  E* 

(3)  The  maxleua  of  A  and  B  Is  defined 

as  AVB  *  {a(T)b:  a  e  A,  bCB  and 
•  -  b}  E. 

(4)  The  aJnloua  of  A  aod  B  Is 


Eight*nelghbor  logic  operations  are  some  of 
the  most  common  operations  used  In  Image  pro* 
cessing*  These  operations  lend  themselves 
particulary  well  to  the  type  of  array  processor 
architecture  portrayed  in  Figure  1*  The  KPP 
accomplishes  eight-neighbor  operations  by  shifting 
over  the  entire  array*  Ve  let  >^(a)  denote  the  3  x 
3  neighborhood  configuration  corresponding  to  this 
wiring,  where  *a"  denotes  the  center  pixel*  The 
links  of  a  processor  to  the  imreedlate  neighbor 
.hoods  can  usually  be  controlled  by  on-off  switch¬ 
es,  allowing  the  configurations  of  different  sub- 
neighborhoods*  It  will  be  convenient  to  label  the 
corresponding  subsets  of  H(a)  by  the  counterclock¬ 
wise  numbering  convention,  illustrated  by  the 
following  f i>  jres: 


M(a) 


Is  defined  as  /.Ih  *  {a(^)b: 
Sc  A,  bcB  and  a  ‘-b|  ^E* 


■*f‘246 


(a)  -  4 


In  contrast  to  tVr  operations  defined  in 
Section  2,  all  the  stvive  operations  arc 
cooDutatlve  and  associative.  In  fact, 

A*0*  -  A,  Axl*  -  A,  A.  A)  -  0*, 


In  part*.cuj*r,  V.  -  7  1®  CLIP'S 
neighborhood  circuitry,  while  ^'■^245  neigh¬ 
borhood  Brrantes*ent  of  the  KPP. 


f 

5 

» 
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II 

(1)  X  -  tS«  *  -  «)  Bud 

T  -  s;  *  “  y} 


(11)  I(u.»)  -  I»p{(-2«l/n)(uX*wY)), 

12 

where,  ii,v  €  Z. 

n 

The  Fourier  neighborhood  function  of 

on 

on  image  la  defined  aa  the  function  FlO 
•  F(S),  where  F^(u,»,0)  -  AjL£(y,v). 

14. 

The  Fourier  cranafonacd  Image,  F(A), 

of  A  coo 

then  be  expressed  by  the  simple  fonsuln 

r(A)  - 

IS 
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